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We investigate the dynamical response of strongly interacting ultra-cold fermionic atoms near 
Feshbach resonance to varying magnetic fields. Following the experimental practices, we calculate 
the response of the atoms to oscillating and to linearly ramped magnetic fields respectively. For 
oscillating magnetic fields, depending on the frequencies and the amplitudes of the oscillations, the 
response of the pair excitation gap shows either linear or rich non-linear behaviour. In addition, both 
the spectral studies through the linear response theory and the time-domain simulations suggest the 
existence of a resonant frequency corresponding to the pair dissociation threshold. For linearly 
ramped magnetic fields, the response of the excitation gap shows damped oscillations. The final 
value of the excitation gap depends on the rate of the field sweep. 



I. INTRODUCTION 

Feshbach resonance in ultra-cold Fermi gases, where 
the energy of two scattering atoms is magnetically tuned 
to coincide with that of a quasi-bound molecular state, 
has provided a rare opportunity to study strongly corre- 
lated many-body systems with tunable interactions 
Recently, considerable experimental efforts have been 
dedicated to the condensation of fermionic atom pairs 
near Feshbach resonance 

U II II II II The properties of 
such a Fermi condensate and its crossover to the Bose- 
Einstein condensation (BEC) of diatomic molecules are 
currently under active studies H, H 0, 0, 0, ^| . 

In this paper, we investigate the dynamical response of 
a Fermi condensate to magnetic field modulations. We 
study two kinds of modulations, the sinusoidal oscillation 
and the linear ramp, both of which have been frequently 
applied in experiments [J 0, 0, H, Q El- There 
have been some recent studies on the atom-molecule os- 
cillations under a sudden switch of the magnetic field 
[Til IT^| . By considering response of the Fermi conden- 
sate to more general variations of the magnetic field, we 
find richer response structures of this strongly interacting 
system under different diving frequencies or sweep rates. 
In particular, oscillating magnetic fields have been used 
recently as important experimental tools to investigate 
the relaxation time of the Fermi condensate [7|; while 
fast magnetic field sweeps from the BCS side of Fesh- 
bach resonance to the BEC side have been used to probe 
the properties of the Fermi pairs on the BCS side 0,3- 

In the case of oscillating magnetic fields, we first de- 
velop a linear response theory in the frequency domain 
to describe the dynamics of the Fermi condensate under 
arbitrary but small variations of the magnetic field. This 
spectral analysis is then compensated by direct time- 
domain simulations. The result shows that the response 
of the Fermi condensate remains linear up to the pair 
dissociation frequency, around and above which the non- 
linearity of the system becomes dominant, giving rich 
non-linear response structure. For linearly ramped mag- 
netic fields, we perform real time evolutions at differ- 
ent ramping rates. The responses of the system are es- 
sentially damped atom-molecule oscillations. The final 



states of the damping process, however, are dependent on 
the ramping rates. For adiabatic ramps, the final states 
are close to the stationary state at the final magnetic field 
of the field sweep. This is consistent with the experimen- 
tal practice, where adiabatic sweeps are applied across 
the Feshbach resonance from the BCS side towards the 
BEC side to efficiently convert Fermions into Feshbach 
molecules jQ, 0|. In contrast, the final state of a sud- 
den ramp correlates rather with the stationary state at 
the initial magnetic field. 

In the following, we first derive the equations of mo- 
tion starting from the general two-channel Hamiltonian. 
Our approach is based on the time-dependent variational 
theory, by assuming an evolving pair wavefunction as the 
ansatz state. Then in Sec. Ill, we apply the equations of 
motion to study the response of the system to oscillat- 
ing magnetic fields. We first perform a spectral analysis 
using the linear response theory to characterize the re- 
sponse spectrum, which is later supported by direct time- 
domain simulations. In Sec. IV, we apply the equations 
of motion to study the response of the system to linearly 
ramped magnetic fields. We conclude with a summary 
in Sec. V. 



II. MODEL HAMILTONIAN AND THE 
EQUATIONS OF MOTION 



The coupling in strongly interacting Fermi gas near 
Feshbach resonance can be described by a two-channel 
model with the Hamiltonian H El E Ellll El : 

H = ^ e k4,a a k, CT + ^(£k/2 + 7)&lA 

k,cr k 

+ ^2 ~ff ( & q a k+q/2,T a -k+q/2,l + h.c) , (1) 
k,q V 

where £k = ?i 2 k 2 /2m is the atomic kinetic energy (m 
is the mass of the atom), g is the atom-molecule cou- 
pling rate, V is the volume of the system, and a£ , 6^ 
are the creation operators for the fermionic atoms and 
the bosonic bare molecules with momentum hk, respec- 
tively. Atoms in two different Zeeman levels are labeled 



2 



by the spin index a( a =f , J.). The bare energy detun- 
ing 7 of the quasi-bound molecular level is related to 
the physical detuning j p = 2/xb [B (t) — B ] by the stan- 
dard renormalization relation 7 = j p + g 2 /V^Z k (l/2£ k ) 
Hi H> 0|- In the expression for y p , B (t) is the applied 
magnetic field which we take to be time-dependent, and 
Bo denotes the Feshbach resonance point. In writing this 
Hamiltonian, we have ignored the background scattering 
as it is negligible near Feshbach resonance. 

The basic method to deal with this system at very 
low temperature is the crossover theory which is based 
on the variational assumption that the ground state of 
the above Hamiltonian can be written as a condensate 
of the dressed molecules [1, EJ 0, ■ Here, we extend 
this variational method to the time-dependent case (with 
varying B (i)) by assuming that the state of the system 
at any time t can still be written in the form: 

|* (<)) = 7VexpL4 (i) bl + Y,h (*) al k ,;4,T] l wac > ' 

k 

(2) 

but with time-varying variational parameters A (t) and 
/k (t) (where J\f is the normalization factor). The param- 
eters A (t) and /k (t) can be identified with (bo) and the 
Cooper-pair wavefunction in k-space, respectively. This 
variational ansatz is reasonable since a condensate state 
of the form of Eq. (2) is pretty robust against excitations 
induced by the magnetic field variations: it is unlikely to 
produce single-atom excitations due to the existence of 
an energy gap (this is particularly the case when the rate 
of the field variation is smaller than the energy gap); it 
is also hard to generate excitations with from the 
q = condensate due to the momentum conservation. 

With the ansatz state above, we minimize the action 
S = Jdt($ (t) I {id t — H) |$ (t))+h.c. under |$ (t)). The 
functional derivatives of S with respect to the variational 
parameters A* (t) and / k (t) give the dynamical evolution 
equations: 

i/k(t) = 2e k f k (t)--?=A*(t)ti(t) + -?=A(t) ) (3) 



iA(t) = ^=V- 



■jA(t). 



(4) 



The above evolution equations can be arranged into 

a more convenient form if one introduces the pseudo- 

— -4- i t v 

spin operators cr k = ak,-fd-k,j,, <?k = -ki a kT' °k = 

7ik,t + ra-k,4. - 1 (where n k ,| = a^ak^HS 113 E3- 
Under the variational state (2), it is easy to check 
that the mean values of these pseudo-spin operators are 
uniquely related to the pair-wavefunction / k (t) through 
/k = ((ok) + 1)/ (2 (<r k )) and the natural constraint 

4( CT k) (°k~) + ( CT k) 2 = !• With these relations, the evo- 
lution equations (3-4) can be written in terms of the ex- 
pectation values of the pseudo-spin operators: 



i(tr k ) = 2e k (a k > - A , 

= 2A( ( r k )*-2A*(a k ) 



(5) 
(G) 



(7) 



where the parameter A = g (bo) = gA (t) is physically 
identified as the excitation gap for fermionic excitations 
[20|. The set of equations (5-7) can also be derived by 
replacing the operators with their mean values in the 
Heisenberg evolution equations for cr k ,cr k , and bo 0, 
Il7| . The derivation here provides a different physical 
interpretation of the basic approximation: it is equivalent 
to an extension of the variational ansatz for the crossover 
theory 0, |j| 0, |2£j to the time-dependent case. 

To investigate how the system responds to the vari- 
ations of the magnetic field (which varies the detuning 
7), we look at the evolution of a detectable property, 
which we choose as the gap function A (t) . The gap A (t) 
is proportional to the occupation in the bare molecular 
level, and is directly measurable, for instance, through 
the radio-frequency spectrosco py | 22l l23| or the bare 
molecule population detection |24j. In the subsequent 
sections, we will examine the response of the system to 
magnetic field oscillations and linear sweeps at different 
rates, respectively. 



III. RESPONSE TO OSCILLATIONS OF THE 
MAGNETIC FIELD 

To characterize the response of the system to arbi- 
trary, yet small, magnetic field oscillations, we expand 
A (t) , (a k (f)) , (tr k (i)) around their stationary values, 
and linearize the resultant evolution equations. The 
stationary values are obtained from the conventional 
crossover theory: we minimize (H — /j,N) with a fixed 
initial detuning 7 P = 70, and get the following gap and 
number equations characterizing the stationary state : 



( 70 -2 M ) = f 7 (^-l---L), 



V^2E k 2e k ' 



2 f + ?(' 



£k - M 



(8) 
(9) 



where N — ^ k g a fc CT ak )(T + 2 6 t &k denotes the to- 
tal atom number, fi is the Lagrange multiplier with 
th e meaning of a tomic chemical potential, and i? k = 
y/ (e k — fi) 2 + Ap denotes the single-atom excitation 
spectrum in the stationary case. From the stationary 
equations, the gap A and the chemical potential [i in 
the stationary state can be solved, from which the sta- 
tionary values of the pseudospins are obtained as follows: 



/ -\ A ° 



£k - fJ- 



(10) 
(11) 



From the linearized evolution equations, we can find 
out the response relation between SA(t) = A(£) — 
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Aoe" 2 ^ 4 and * 7 (t) = 7p - 7o = 2fi B [B (t) - B (0)], 
where B (0) is the initial magnetic field. It is conve- 
nient to express this response relation in terms of their 
Fourier transforms 5A(ui) = J SA^e^'^dt and 

I 



Sj(uj) = J Sj(t)e~ luJt dt. The general response relation 
can then be written as 8A(u>) = P (u>) 6j(u)), with the 
response function P(u>) given by: 



p( x Ap (FH-G*(-o;)+a;-(7 -2 At )) 

1 > ( W + G(w) + ( 7o -2/i))(- W + Gf(-w) + ( 7o -2^))- J FHF*(-w)' 1 J 



in which 



2w 2 E F J 4E k (El-f%y v ' 2^ 2 £W v AE k {E 2 k - £) 2e fc 



0.02 




FIG. 1: (Color online) The response function shown in the 
range \u\ < 2Ao (here, 2Ao ~ 0.87Ef, and u> is in the unit 
of Ef)- Other parameters are comparable with the 40 K ex- 
periments 0, where gy^ ~ 13.8Ef (n is the atom number 
density and 70 ~ 100-Ef). 



The response function P(lu) specified by the above 
equations can not be directly compared with the experi- 
mental measurements as it is |A| that is measured. It is 
easy, however, to construct the response function R (u>) 
for S\A\ = \A\ — Aq from P (u) with the following simple 
connection: 



R{u) = -[P{w)+P*{-w)]. 



(13) 



The resultant spectrum for R (to) is numerically calcu- 
lated and shown in Fig. 1, with u> in the range — 2Ao < 
lo < 2 Ao , and other parameters comparable with the 
40 K experiments |3J. The spectrum shows that the re- 
sponse of the Fermi condensate is linear for small driving 
fields with the response function only weakly dependent 
on the driving frequency. The spectrum also shows res- 
onance behaviour at u) = 2Aq, which can be associated 
with the pair dissociation frequency. 

To complement the above linear response calculations, 
we also simulate the evolutions of the gap function A (t) 
from Eqs. (5-7) directly in the time domain, with si- 



nusoidal magnetic fields at different driving frequencies. 
With this temporal simulation, we would like to first 
confirm the resonance predicted by the linear response 
theory and second find out the validity region for the 
linear response theory and reveal some non-linear be- 
haviours. In the simulation, we start from an equilib- 
rium state at the positive detuning (B (0) > Bo), and 
assume that the variation of the detuning takes the form 
8j(t) = 7p - 7 o = 7am P sin (uot), with 70 = 100£> (corre- 
sponding to (B (0) - B ) - 0.2G for 40 K) and a small os- 



cillation amplitude j amp — lEp (Ep = (S7T 2 h 



2*3 x2/3 



/2m 



is the Fermi energy, where n is the atom number den- 
sity and m is the mass of the atom). The results of the 
evolutions are shown in Fig. 2 with several characteristic 
driving frequencies. For driving frequencies much smaller 
than 2Ao (Fig. 2a), the <5|A(i)| follows the same sinu- 
soidal curve as (B (0) — B (t)) with almost no detectable 
delay. This shows that the system responds very fast to 
the changes of the magnetic field, with a time scale much 
shorter than the thermal relaxation time reported in the 
recent experiment 0, l2o| . When the driving frequency 
approaches 2Ao (Fig. 2b), some non-linearities show up 
in the response with a periodic modulation of the oscil- 
lation amplitude. For u> — 2Ao (Fig. 2c), the resonance 
behaviour is apparent. The amplitude of the oscillation 
continuously increases within the time interval attainable 
in the simulation. For the driving frequency significantly 
larger than 2Ao (Fig. 2d), the response is almost linear 
at the beginning, then the nonlinearity shows up causing 
increase of the oscillation amplitude. We have also sim- 
ulated the time evolution of the gap function A (t) with 
a much larger oscillation amplitude of the driving field 
corresponding to 7 amp = 50Ep. For small non-resonant 
driving frequencies, the response is still linear even with 
such a large driving amplitude. For instance, with the 
same u> as in Fig. 2a but 7 amp = 50Ep, the response 
curve of 6\A(t) | still has the same shape, although its 
amplitude is amplified by a factor of 50 as one expects 
from the linear response theory. For driving frequencies 
larger than 2|Aq|, even if they are non-resonant, rich 
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FIG. 2: (Color online) The real-time (in the unit of 1/J5f,~ 120/is for 40 K experiments Q) evolutions of the pair excitation 
gap |A(t)| (in the unit of Ef) shifted by the initial gap Ao (Ao ~ 0.433-Ef )• The system is driven by sinusoidally oscillating 
magnetic fields at different driving frequencies to. The solid lines represent the evolutions of the excitation gap 5[A(t)[, while 
the dashed line is the re-scaled magnetic field variation B(0) — B(t). In all cases, ^ amv = IEf, Qsfn ~ 13.8Ef, and 70 = 100-Ef- 
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FIG. 3: (Color online) The real-time evolution of the pair 
excitation gap |A(t)| shifted by the initial gap Ao. The same 
parameters as in Fig. 3d but with a larger driving amplitude 
Tamp = 50-Ef, which roughly corresponds to 0.15G in the 40 K 
experiments Q]. 



non-linear structure shows up in the response, as one 
can see in Fig. 3. 



IV. RESPONSE TO LINEARLY RAMPED 
MAGNETIC FIELDS 

Another branch of important experimental practice is 
to probe the properties of the Fermi condensate by ramp- 



ing the magnetic field from the BCS side of the Feshbach 
resonance to the BEC side. Slow adiabatic ramps are 
used to create molecular condensates 0, , while fast 
sudden ramps are used to project the correlated pairs 
onto Feshbach molecules so as to probe the properties 
of the Fermi condensate on the BCS side of the reso- 
nance It is desirable then, to understand the ramp- 
ing process and the effects of the ramping rate. In this 
section, we study the response of the system to linearly 
ramped magnetic fields with different ramping rates. 

Similar to the previous section, where we vary the mag- 
netic field by modulating the detuning in Eq. (7), now we 
impose a detuning in the form 7 P =70— (jo — Jf)t/T ramp , 
where 70,7/ are the initial and final detunings respec- 
tively, and T ramp is the duration of the field sweep. To 
facilitate comparison between processes with different 
ramping rates, once the detuning has been ramped to 
7y, the final detuning, we let the system evolve under 
this constant magnetic field. Again, we use the pair ex- 
citation gap A(t) to represent the response of the Fermi 
condensate. As the numerical simulation is quite time 
consuming at large coupling rates, for the sake of ef- 
ficiency, we set the atom-molecule coupling rate to be 
g\fn ~ IQEf, which is slightly smaller than that of the 
40 K system, where gy/n ~ 13.8Ep3\- Qualitatively, the 
results should not be affected by this change of coupling 
rate. We perform real-time evolutions of the excitation 
gap based on Eqs. (5-7), and the results are shown in 
Fig. 4. The most prominent feature of the results is 
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FIG. 4: (Color online) The real-time evolutions of the pair 
excitation gap |A(f)| under linear magnetic held sweeps with 
different sweep rates. While the thin solid curve repre- 
sents the evolution under a sudden ramp, the ramping time 
T ram p, as defined in the text, for the rest of the curves are: 
Tramp = 1 / Ef for the thin dashed curve, T ra mp = 5/Ep 
for the thin dash-dot curve, T ra mp = W/Ef for the thick 
solid curve, T ra m P = 15/ Ef for the thick dashed curve, and 
Tramp = 20/ Ef for the thick dash-dot curve. The atom- 
molecule coupling rate is taken to be g\/n ~ WEp, with 
the initial detuning 70 = WOEf, and the final detuning 
7/ = -500E F . 
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FIG. 5: (Color online) Comparison between the final excita- 
tion gap values resulting from slow ramps (T ra m P = 20/ Ef) 
to different final detunings, with those of the stationary states 
at the corresponding final detunings. The solid line represents 
the stationary gap values, and the dashed line represents the 
final gap values of the slow ramps. The starting points of all 
the field ramps are fixed at 70 = IOOEf, while the final de- 
tunings are varied over the range of 7/ = ~ — 500-Ef- The 
atom-molecule coupling is taken to be g^/n ~ WEf- 



that the excitation gap undergoes damped oscillations 
under fast field sweeps. This result is in agreement with 
Ref. for the particular case of a sudden sweep, but 
is in contrast to the undamped oscillations reported in 
[l2|. We expect that this kind of damped oscillations 
is a general response behaviour for a system with many 
degrees of freedom as is the case here; the undamped os- 
cillations on the other hand should be a more artificial 
effect resulting from some unrealistic constraints. For 
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FIG. 6: (Color online) Comparison between the the final ex- 
citation gap values resulting from sudden ramps from differ- 
ent initial detunings, with those in the stationary states at 
the corresponding initial detunings. The solid line represents 
the stationary gap values, and the dashed line corresponds 
to the final gap values of the sudden ramps. The final de- 
tunings of all the field ramps are fixed at 7/ = — 500-Ef, 
while the initial detunings are varied over the range of 70 = 
100-Ef ~ — WOEf- The atom-molecule coupling is taken to 
be g\pn ~ 10-Ef- 

slower sweeps, the damping also occurs along with the 
sweep but with a relatively smaller amplitude; and if the 
sweep rate is slow enough, the absolute value of the exci- 
tation gap will follow the field sweep almost adiabatically. 
Another important feature is that the final state of the 
system, reflected in the absolute value of the final gap, 
depends on the rate of the field sweep. Even though the 
initial and the final detunings are set to be the same in 
all the evolutions in Fig. 4, there appear to be consid- 
erable differences in the values of the final gap at dif- 
ferent ramping rates. These rate-dependent "stationary 
states" are different from the stationary states calculated 
in Sec. Ill with Eqs. (8-9). The final state at a given 
ramping rate is one such that the absolute value of the 
excitation gap is not evolving but its phase changes with 
time. Even under such a final "stationary" state with a 
constant |A (t)\, numerical analysis shows that there still 
appear to be complicated momentum-dependent oscilla- 
tions of the pseudospins (c^), (ck)> an d (°k)) which 
effectively describe the shape oscillations of the dressed 
molecules. The amplitude of such oscillations decreases 
as the sweep rate approaches the adiabatic limit, and in 
the adiabatic limit, the state just follows the real station- 
ary state of the system determined by Eqs. (8-9) as the 
detuning decreases. 

Among these many ramping-rate dependent final 
states, of particular interest to us are the final state of 
a sudden ramp and the final state of an adiabatic ramp, 
due to their close experimental relevance. To study the 
property of the final states, we compare the excitation 
gap in the stationary state at either the initial or the final 
detuning with the gap in the final state of the evolution. 
One can either fix the initial detuning while changing the 
final detuning of the field sweep or the other way around. 
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We have plotted in Fig. 5 and Fig. 6. the results of our 
simulations. For a very slow ramp, if one fixes the ini- 
tial detuning and varies the final detuning, the final gap 
would approach the stationary values of the gap at those 
different 7/. The differences between the two curves in 
Fig. 5 can be made smaller by decreasing the ramping 
rate further. This is in agreement with the experiment 
|14|. where slow field sweeps, which roughly correspond 
to ramp time of over 100/ Ep to cover the same detuning 
range in Fig. 4, were applied to create molecules. In a 
sudden ramp, one needs to fix the final detuning while 
varying the initial starting point. It appears, as in Fig. 
6, that there exist correlations between the final gap and 
the initial gap, while both are much smaller than the 
stationary value of the gap at the final detuning. There 
are no such correlations if one fixes the initial detuning 
and vary the final detuning. This result is in agreement 
with Ref. from a single channel calculation. The 

finding also agrees with the experimental practice that 
a fast projection would allow us to probe information in 
the initial state. 



V. SUMMARY 

In summary, we have investigated the dynamical re- 
sponse of a Fermi condensate near Feshbach resonance 
to external magnetic field modulations. For a sinusoidal 
modulation with small amplitude and frequency, the re- 
sponse of the Fermi condensate linearly follows the mod- 



ulation. The linear behaviour breaks down when the fre- 
quency of the modulated field approaches the pair disso- 
ciation frequency 2Ao- The response of the system will 
be composed of highly non-linear structures if the fre- 
quency increases above this threshold. Calculations of 
the response spectrum using the linear response theory 
as well as real time evolutions consistently support the 
conclusions. We have also calculated the real time evo- 
lutions of the pair excitation gap under linearly ramped 
magnetic fields, the results of which are consistent with 
the experimental observations. The general feature of the 
response to the field sweep is the damped oscillations, 
which are more manifest in faster sweeps. The final state 
of the sweep depends on the rate of the sweep. While the 
final state of a slow ramp is quite close to the stationary 
state at the final detuning; that of a sudden ramp corre- 
lates with the stationary state at the initial detuning, as 
can be seen from the correlations between the final gap 
and the initial gap in a sudden ramp. The results are in 
agreement with the experimental practices that adiabat- 
ically slow field sweeps efficiently convert Fermion pairs 
on the BCS side of the Feshbach resonance into bound 
molecules deep in the BEC side; and that sudden field 
sweeps project the Fermion pairs onto the molecules, so 
that the properties of the initial state can be probed via 
measurements of the final state. 
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